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We study a readout scheme of superconducting flux qubit state with a Cooper pair box as a 
transmon. The qubit states consist of the superpositions of two degenerate states where the charge 
and phase degrees of freedom are entangled. Owing to the robustness of transmon against external 
fluctuations, our readout scheme enables the quantum non-demolition and single-shot measurement 
of flux qubit states. The qubit state readout can be performed by using the non-linear Josephson 
amplifiers after a 7r/2-rotation driven by an ac-electric field. 
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I. INTRODUCTION 

High-fidelity detection schemes has been intensively 
studied to reduce the decoherence during the readout 
process of qubit state. The dispersive measurements- 
is known to minimally excite the spurious degrees of 
freedom from environment and has low backaction on 
the qubit. For superconducting qubits this measure- 
ment has been performed for circuit quantum electro- 
dynamics (QED) architecture 4 - - — , quantronium qubit 7 ", 
and superconducting flux qubits*. In order to detect 
the qubit states, however, the readout process should 
be fast compared to the qubit relaxation time and not 
invoke the transition between qubit states. For the 
transmon qubit^r— this kind of dispersive readout has 
been implemented by using bistable hysteretic system 
of non-linear resonator such as the Josephson bifurca- 
tion amplifier— and Josephson parametric amplifier 13 . 
The readout by the non-linear Josephson resonator en- 
ables the single-shot readou t 12 ' 13 and the quantum non- 
demolition (QND) measurementi 2 - - — for transmon qubit. 
The measurement of qubit states in a single readout pulse 
is mostly important for the scalable design of quantum 
computing. For the error correction in quantum algo- 
rithm code the fast and efficient single-shot readout is 
indispcnsible. Moreover, the transmon qubit remains 
robust against relaxation during the measurement and, 
thus, the eigenstate of qubit is not changed, which en- 
ables the QND measurement. 

In this paper, we propose a new scheme for readout of 
flux qubit states coupled with a Cooper pair box. The 
qubit states consist of the current states of the flux qubit 
loop, and are manipulated by a magnetic microwave. The 
phase degree of freedom of flux qubit 1 - 5 -"— and the charge 
degree of freedom of the Cooper pair boxi^r— are entan- 
gled with each other so that the flux qubit state may be 
read out by detecting the charge state. Since the phase 
and charge variables are canonically conjugate with each 
other, these variables cannot be determined simultane- 
ously After rotating the qubit state with an oscillating 
electric field at the end of the qubit operation, the qubit 
state measurement can be achieved by charge detection. 

Figure [1] shows our design for readout of flux qubit 



states. The upper part is the three-Joscphson-junctions 
qubit (flux qubit) and the lower part consists of a large 
Josephson junction (-E74) and a Cooper pair box between 
two small Josephson junction (Ej 5 and Ej 6 ). The lower 
part is similar to the quantronium qubi t 24 ' 25 . For the 
quantronium qubit the charge state is detected by mea- 
suring the output pulse depending on the phase difference 
across the large Josephson junction. On the contrary, our 
qubit design aims to read out the flux qubit state of the 
upper part by detecting the charge state in the lower 
part. 

In the present design we consider the transmon qubit 
with a large shunted capacitance C s as a Cooper pair box. 
Owing to the large shunted capacitance the transmon is 
robust against the charge fluctuation. The qubit state 
measurement can be performed in a dispersive manner by 
using the Josephson non-linear resonators. By detecting 
the transmon state we will be able to read out the flux 
qubit state in a non-destructive single-shot measurement. 
On the other hand, the optimal point measurement of 
flux qubit states has been studied previousl y 26 ' 27 . In our 
design also the qubit state readout can be performed at 
an optimal point. 



II. HAMILTONIAN OF COUPLED QUBITS 

The total energy of the system in Fig. [T] consists of 
the Josephson junction energy and the charging energy, 
neglecting small inductive energy. The Josephson junc- 
tion energy Ujj({4>i}) is represented in terms of the phase 
differences (f>i across the Josephson junctions, 



U 



j j 



({<&}) = ~^2 E Ji cos < 



(1) 



We have two boundary conditions for the upper and lower 
loops, 



</>2 + </>3 + 04 = 27r(n a + /), 
04 + <fc + </>6 = 2Tf(n b - f h ) 



(2) 
(3) 



with integers n a and n^. Here, / = ^cxt.a/'I'o an d 
fb = ^cxt.b/^o with the external fluxes <E> xt,a and $ e xt,b 
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threading the upper and lower loop, respectively, and 
= h/2e. We, for simplicity, set 



Ej2 — Ej3 

and thus we have 



Ej, Ej 5 — E, 
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E 



lb, 



2irm a , 
2i:m b 



(4) 



(5) 
(6) 



with integers m a and m . Hereafter, we set Eji — 0.8Ej. 

Since we have 4 constraints of Eqs. ©-(JU), the Joseph- 
son junction energy in Eq. ([1} can be represented in the 
plane of (0 a , D ) as shown in Fig. [2J where 



0a = (0 2 +03)(mod27r)/2, 
4>b = (05 + 6 )(mod 2tt)/2. 



(7) 

(8) 



In Fig. H we set / = 0.5 and f b = 0.5, and | 1) (| t>) de- 
notes the counter-clockwise (clockwise) current state in 
the upper loop. We can observe that the current direc- 
tions of upper and lower loops are correlated with each 
other. When a > O(0 a < 0), b > 0(<j) b < 0), which 
means, depending on the current direction in the upper 
loop, the sign of phase shift across the junctions 5, 6 
changes. Then we can set 



; (mod 2tt) = 6 (mod 2tt) = ±0/2 



(9) 



for | 1) (| f)) state with = 2|0 b |. 

On the other hand, the transmon qubit states are de- 
scribed by the number of Copper pairs, |n) and \n+ 1), 
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FIG. 1: Flux qubit with a large Josephson junction of Joseph- 
son coupling energy Ej4. A Cooper pair box with a large 
shunted capacitance C s (transmon) is attached to the flux 
qubit across the large junction. 
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FIG. 2: The Josephson junction energy Ujj({(f>i}) in the 
plane of (4> a ,(j>b), where cf> a — (4>2 + <^>3)(mod 2n)/2 and 
06 = (05 + 0e)(niod 2n)/2. At the local minima of Ujj, cj> a 
and <f>b have the same sign. For the flux state | 4-), the per- 
sistent current in the upper (lower) loop is counter-clockwise 
(clockwise) as shown in Fig. [T] and vice- versa for | f). 



in the Cooper pair box (denoted as B in Fig. [T]). 
charging energy splitting^ is given by 



E., 



E n 



2E ch = -4E c (l-2n g ) 



The 



(10) 



with the charging energy Ec = (2e) 2 /2Cs of the trans- 
mon and the dimensionless gate voltage n g = C g V g /2e. 
Here, V g is the gate voltage and Ce = C g + Cj + C s with 
the gate capacitance C g and the Josephson junction ca- 
pacitance Cj for transmon. The transition between the 
states | re) and |re+l) is invoked by the Josephson junction 
energy Ej^^.. 

The Hamiltonian Hq for our qubit with / 0.5 and 
f b w 0.5 is given by 

H = I®[-E ch \B)(B\ + E ch (\L){L\ + \R){R\)} (11) 
+(-e| I) (I | +e| t)(t I - AUXt | - A| t>(4 |) ® / 

-I I ® (e^\R)(B\+e^\B){L\ + c.c. 
-I t)(t I ® ^ (e-~^\R){B\ + e-**\B){L\ + c.c.) 

The first term shows the energy levels of the Cooper pair 
box. When an additional Cooper pair is in the box, the 
state is represented as \B), while if the Cooper pair tun- 
nels into left or right side of the box, the state is denoted 
as \L) or \R). For the state \B), the number of Cooper 
pairs in the box is re + 1, while for the states \L) and \R) 
it is n. The second term describes the dynamics of the 
flux qubit, where A is the transition rate between | 1) 
and | f), and e = <J>o/ p <5/ is the energy level shift of flux 
qubit part with 5f = f — 0.5 and I p being the persistent 
current in the flux qubit loop. The latter two terms de- 
scribe the tunneling of a Cooper pair. The phase shifts 
involved in the tunneling have different signs depending 
on the flux qubit states. 
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We consider the case that Ej4 ^> Eji similarly to the 
quantronium qubiliSi. As a result, 04 is very small com- 
pared to other (pi's, which allows an analytic analysis 
for Hq. Then, from the boundary condition of Eq. ([3]) 
around the lower loop, we obtain <j> rj 7r for f b = 0.5. 
If we set (j) = tt, we can obtain the eigenvalues of the 
Hamiltonian Hq as 



Ao,i,2, 3 = ±\/F±2VG, 

^4,5 



±y/A 2 + e 2 /2 



(12) 
(13) 



^ + A 2 + E 2 h +Ej b /2 and G = e 2 E* h + £ c 2 h A 2 - 



with F = 
e*E%/2. 

For an analytic analysis we set e = (/ = 0.5). Then 
the eigenstates becomes 



IVo) = N a 



|Vi 



|V>3 



iVi 



AT. 
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(U) - 1 t»(|i) 
~— (-Q 

E Jb vv ^" 

(U> + lt»(|i) 
— (-Q 

E Jb vvr 

(U> + lt»(|£) 

2i ( 1 



£jb \V2 



-A+J (U) + |t»|B) 

(14) 

-A_) (\i)-\t))\B) 

(15) 

t»|S> 
(16) 

t»|B) 
(17) 



+ A_) (||) 
I*))]. 

" (II) 



+(|;)_|t»(|i>-|i2»]. 
with normalization factors JVj and eigenvalues 

Ao = --j=Q + ,M = --^Q-, 

A 2 = —Q_,X 3 = —Q + , 

and 



<9± = y E jb 

A± = A±LE ch 



2A 2 fc) 



The excited states ^4), jV's) corresponding to A4 
and A5 = A + can also be obtained as 

m = -j=(\u + \t))(\L) + \R)), 

|V>5> = ^=(|1>- 1 1>)(|£> + !#>)• 



(18) 
(19) 

(20) 
(21) 

-A_ 

(22) 
(23) 



In the eigenstates \ipo), \ipi)i \1fa1) an d | ^3} , the flux states, 
I I) + I f) and I 4.) — I f), and the charge state, \B) and 
\L) — \R), are entangled with each other, whereas the 
states, \%l>4) and {ips), are product state. 

In Fig. [3(a) we plot four eigenvalues Ao, Ai, A2, A3 
in the plane of (f,n g ). Here we use Ej/Eq = 100 for 
transmon. Figs. [3(b) and (c) show the cut view of Fig. 



[3(a) along / with n g = 0.5 and along n g with / = 0.5, 
respectively. Here, the point (f,n g ) = (0.5,0.5) is an 
extreme point. At this point the states are degenerate 
and we will define the superposition of these degenerate 
states as a qubit state later. 

The Rabi oscillation between the eigenstates can 
be performed by applying a magnetic microwave field 
£u)(t) = g cos uit on the qubit, where g is the coupling 
strength between the qubit and the microwave. The to- 
tal Hamiltonian is given by H = Hq + H mw , where H mw 
describes the interaction between the qubit and the mi- 
crowave, 

H mw = (-e u (t)\ I) (I I + <*,(*) I t)(t (24) 

We transform the total Hamiltonian H to H in the basis 
{\if) ), IV'i), 1^2), 1^3), |-04), |-05)} as, 



H= 



— Si cos 52 cos cut 
-gi coscdt — -^=Q_ — <7 2 cosojt0 
g2C0suit "72^" —gicosuito 
— c/2 cos uit —gi cos wt 










\ 








A_ 

A + y 



(25) 



with 



9i 



!J2 



2V2A + Q- + Q+ Q--V2A- 



2y / 2A - Q_ 



q+ + \/2A H 

lQ- + V2A, 
Q+ + V2A^ 



.9, (26) 



■9- (27) 



Energy (E.) 
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FIG. 3: (a) Qubit energy levels, Arj,Ai,A2 and A3 with 
A/Ej = 0.001, Ejb/Ej = 0.1 and E J4 /Ej = 50. The Rabi 
oscillation is performed between the local extreme points of 
energy planes. Cut views of (a) with (b) n g = 0.5 and (c) 
/ = 0.5. 
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Since the states, \i(>a), IV's), are decoupled from the 
other states in H, we hereafter truncate these states 
from the basis {\tpo), [fa), IV^J^i), l^s)} and 
represent the Hamiltonian matrix H in the basis 

{\iM,\4>i),\ih),\ik)}- 

In order to obtain the transition probability, we intro- 
duce a rotating frame such that cf>(t) = T)ip(t). Then 
the Schrodinger equation Hip(t) = ih-^tp(t) is written as 
if t (j){t) = H R 4>{t) with 



which shows a Rabi oscillation between |^o) 
with the Rabi frequency 



and 



M 

UJQ 



E 2 



A 2 . 



and |^ 3 ) 



(35) 



(36) 



H R = VHV^ - W—VK 
dt 



(28) 



We can choose the transition matrix T> between I?! and 

£>2, 



£>i= 
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—iujt 








^ 


















1 












„ — icut 














1 


i 

































1 



(29) 



By choosing T>i(T>2), we can calculate the transition 
probability between the states \i/jo) and \tps} (| Vx) and 
\il>2)). The Hamiltonian H R for the transition between 
|V>o) and ^3) becomes 



H R = 



-f(l 
f(l 



-U) 
J2iujt ^ 



-f(i 



f(i 



71 ( 









2iuit^ 
-f(l 





— 2?u;i 



(30) 



V 







f (1 



e 



f(l 



2 



73* 



Here, note that the Rabi oscillation can be performed by 
the resonant microwave with the frequency 



1 



(31) 



'+• 



not 2 x 

At the operating point (/, n g ) = (0.5, 0.5) the Rabi os- 
cillation between |^»o) and 1^3) may be analyzed in the 
rotating wave approximation (RWA). The usual RWA 
neglects fast oscillating mode, so we set e ±2%tJJt = 1 in 
the Hamiltonian H R —. At this point, Q + = thus 
the eigenvalues arc degenerated and the transition am- 
plitudes are reduced as 



.91 



.92 = 



V2A 



Ejb 



E 2 



E 2 ' 

^Jb 



(32) 
(33) 



In this case the Rabi oscillation demonstrates the max- 
imum fidelity (F=l). If one of g\ and (72 vanishes, we 
were able to check that the fidelity is zero, which means 
that the Rabi oscillation involves two stages. 

In the Hamiltonian of Eq. Q30p there are two ways in 
which the initial state |?/>o) evolves to the final state ^3). 
Note that there is no direct transition amplitude between 
|V'o) and I ^3) states. First of all, the initial state can 
evolve to the final state through the intermediate level 
with energy -(1/V2)Q-. The transitions can be done 
by the off-diagonal terms as 



IV'o) -> IV'l) -> |^3>, 



(37) 



where the first and second steps are driven by the terms 



) and -f (1 



±2iu>t 



) in the Hamiltonian 



-f(l 

of Eq. (|30|) . respectively. The second way is to pass by 
the level with energy (l/-\/2)Q- such that 



H'a) -> U'2) -> 



(38) 



where the first and second steps are driven by the terms 
f (l + e ±2lw *) and -f (1 + e ±2,; "*), respectively. As a 
consequence, the Rabi oscillation is performed in a two 
way double stage manner. It can easily be checked that 
two matrices H R describing these processes in terms of 
either T>\ or T>2 commute with each other, so the total 
unitary evolution is represented as the product of these 
two evolution matrices. 

As shown in the above evolution e~ tH t <p(0) of Eq. 
(|3"4"|) . even though the energy levels of the states |V>2) and 
1^3) are degenerate for n g — 0.5, the initial state |^o) 
always evolves to \tp3), not \ip2)- Numerically also we 
were able to confirm the deterministic relation, 



\ipi) -> \ipa) 



(39) 



through T>\ and T>2, respectively. This means that the 
Hilbert spaces spanned by the basis {\ipo}, ^3)} and by 
the basis {l^i), IV^)} are effectively decoupled. 

In this case the qubit states, |£o) and are the 

superpositions of these degenerate states, 

|£o) =o# )+/^i), |6>=a|^>+^3>. (40) 



We found that for 31 = 32, i.e., Ejb = %/2A, and g <C u>, 
the initial state </>(0) = (1,0, 0, 0) evolves such that 



-^*V(0) = (cos M tj 0, 0, sin M t ), 



ui 



(34) 



The values of a and (3 are determined when the qubit 
states are prepared. A Rabi pulse generates a superposi- 
tion between the qubit state, 



l*(0)) = j>|&) + «lfi 



(41) 



5 



We checked numerically that the single qubit phase evo- 
lution can be achieved by the Larmour precession, 



\^(t))=e ie p\^) + e- ie q\^). 



(42) 



The single qubit Rabi oscillation can also be demon- 
strated. In Fig. |4] (a) we plot the fidelity of the Rabi 
oscillation with the initial state 1^(0)) = |£o), where the 
fidelity is defined as the overlap F = \(^(fl R t = 7t)|£;l)| 2 . 
Here we fix the resonant frequency u) = Q + / \/2 = ujq 
as in Eq. (|31[) . The Rabi oscillation frequency is nu- 
merically obtained as fl R /2n m 70MHz for A/Ej=0M 
and g/Ej=0.05 with Ej/h=100GKz. At the operation 
point (/, rig) = (0.5, 0.5), the fidelity error SF = 1 - F ~ 
O(10" 3 ) for A/£;j=0.005, which is sufficiently small. 

Figures 0] (b) and (c) are cut view of Fig. 0] (a). 
In Fig. 0] (b) the peak width is much broad because 
the Cooper pair box is the transmon, while the peak 
width at F = 0.99 in Fig. 0] (c) is Sf ~ 5 x 10~ 5 
for A/i?j=0.002. Flux fluctuation is estimated to 
be ~ 10- 6 [$ /Hz 1/2 }^i2i and flux amplitude can be 
controlled up to the accuracy of 10 _5 <I>o. Hence both 
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FIG. 4: (a) Fidelity of Rabi oscillation with A/Ej = 0.002, 
Ejb/Ej = 0.1 and Ej^/Ej = 50. Cut views of fidelity with 
(b) / = 0.5 and (c) n g = 0.5. 



the peak widths are sufficiently large. The fidelities 
are higher for smaller value of A/Ejb in numerical 
calculation, whereas F has maximum at Ejt, = \/2A 
in the RWA. The reason for this discrepancy is that 
the oscillating term as well as the constant term 
in (1 + e ±2lwt ) in H R contributes to the transition 
through the intermediate level. Even for the param- 
eter regime where the RWA works well such that 
g/Ej = 0.001 < u>, the fidelities F RW a from the RWA 
and -F n from the numerical calculation arc different 
from each other; we have, for example, (-Frwa,-Fii) = 
(0.735, 0.848), (0.889, 0.782), and (0.973,0.713) for 
A/Ej = 0.04,0.05, and 0.06, respectively, with 
Em/Ej = 0.1. 



III. QUBIT STATE READOUT 

The flux qubit state has usually been detected by 
using the superconducting quantum interference device 
(SQUID). The SQUID with resistivcly shunted Joseph- 
son junction is biased by a current a little below the crit- 
ical current so that depending on the flux qubit state the 
SQUID turns into voltage state to produce the readout 
result, giving rese to decoherence due to the backaction 
to the qubit state. The dispersive readout scheme with 
the non-linear Josephson resonator reduces the backac- 
tion and thus enables the non-destructive measurement. 
For the transmon qubit the single-shot readout has been 
performed through the dispersive measurement schemes. 

For our qubit, which is a hybrid of the flux qubit and 
the transmon, the flux state can be read out by detect- 
ing the state of transmon. Hence, in the present scheme 
the non-destructive single-shot readout can be achieved 
by using the dispersive measurement schemes such as the 
Josephson bifurcation amplifier and the Josephson para- 
metric amplifier. Owing to the large shunted capacitance 
the transmon is robust against the charge fluctuation, 
whereas the qubit operation time becomes long due to 
the flat energy band. Since in our study the transmon is 
not used as a qubit, rather as an auxiliary readout ele- 
ment, the qubit operation time is determined by the char- 
acteristics of the flux qubit part. Hence, our qubit read- 
out design has the advantage of increasing the transmon 
capacitance sufficiently, not worrying about long qubit- 
opcration time. 

We can observe in Eqs. (fl"4"|) - (|17p that the qubit states 
l^i) do not have definite magnetic moment or charge 
number. At the operating point n g = 0.5, l^o) and | "03 ) 
in Eqs. (fT4")) and (jTTJ) approaches the states, 

!*,> - V l) -J t} lL) -J R) + -j= U) + ' t) |g),(43) 



s/2 V2 sft V2 V2 

1 \j)-\t)\L)-\R) i U) + |t) 
y/2 \/2 V2 y/2 V2 



B) 



as A/Ejb decreases. For both states the probability for 
having an additional Cooper pair is 0.5, which means 
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we cannot discriminate between \ipo) and ^3) through 
charge detection. 

In order to read out the qubit state, we apply an os- 
cillating electric field on the qubit, 

H osc =I®g e cosco e t[-(\L)(L\ + \R)(R\) + \B)(B\].(U) 

In the basis {|V>o}, \i>2), \ip3)}, the transformed 
Hamiltonian H osc is represented as 



g e cos u> e t 



( V2A + 

Q+ 
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(45) 



Since Hq in this basis is diagonal, the total Hamiltonian 
Hrc&d = Hq + H osc is transformed to a direct sum of 2 x 2 
matrices, 



-Bread = H Ti 



ead ^ n read ' 



(46) 



where the basis is {|V>o), ^3)} for #r ea d and { I V^x) , ^2)} 
for £T 2 cad . By a proper rotation in the 2x2 subspace the 
qubit state can be transformed to a state which have a 
definite charge number. 

For an analytic analysis we consider n g = 0.5 (E c h = 0) 
and negligible A/Ejb, then 



u 1 

-"read 



TT 2 

n read 



//j 



"7f 



°2 — 9e COS U) e t(J x 



(47) 



The transition between {|^o)) l^)} is described by 



In a rotating frame the Hamiltonian is given by 

Eq. ((281) with X> = e iff * <"<=*. Further, in the RWA the fast 
oscillating modes are neglected as e ±IWet = 1, resulting 
in 



read ' 
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read 
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If we set w e resonant with the qubit energy gap 



2 x E. 



lb 



(48) 



(49) 



by using the time evolution of the qubit state \ip{t)) = 
e~^rcad*| ? /i(0)) the evolutions of initial states \iPo{0))r 
and \ip3(0))n in the rotating frame after 7r/2-rotation 
{geT r = 7r/2) are given by 



|^(0))« = 



\M*/2))R = -j= 



(50) 
(51) 



If the condition uj e T r = 2tt(p — 1/4) with integer p 
is satisfied, the final states in laboratory frame, \ip) = 
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FIG. 5: Readout fidelity F of the state ^ (tt/2)) after tt/2- 
rotation. Here, Ejb/Ej = 0.1 and Eja/Ej = 50. 

e lcr =" c '|^)ij, have definitely different charge numbers as 
follows, 



!*(§)> 



1 



l(|^o(0)) - |^ 3 (0)>) 



^(iv^o(o)> + mm 

l + i|4.)-|t) \L)-\R) 



(52) 



(53) 



s/2 V2 V2 ' 

where |Vd(0)}, |Va(0)) are |^ ), |Va) in Eq. O, re- 
spectively. This condition with p = 10 and g e T r = tt/2 
implies that the coupling constant should be adjusted 
as g e ps 360MHz with £ J6 = O.lEj and £;, 7 =100GHz. 
We also carry out the analysis for the transition between 
{1^1)1 IV'a)} described by -ff 2 cad , resulting in 



l^i(f)) 



\M~)) 



1 



1(1^(0)) - |^(0)» 



V2 
1 + i 



B) 



1 + i I I) -It) 

(|Vi(o)) + |V 2 (o))) 
i + i|4.)+|t) |£)-|ii> 



(54) 



(55) 



V2 \/2 
If we change the basis into the Cooper pair number state, 

\L) \R) 



\n + l) = \B), \n) 



V2 



(56) 



the qubit states of Eq. (140)) are written as 
l£o(f)> = [^(U) + |t» + ^=(U>-|t))] \n + l), 
lag)) = [-^(| I) + 1 1» + 4(1 I) - I t))j In). (57) 



In this case the readout fidelity is defined as the overlap 
with the number state, 



F 



l(£o(f)|n + l)| 2 or F = \(^(^)\n)\ 2 . (58) 
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In Fig. [5] we show the numerical results for the fidelity 
F of |£ (tt/2)) state, where the fidelity error SF = 1 — F 
is as much small as SF ~ O(10" 4 ) for A/Ej = 0.001. 
Non-destructive single-shot charge detection can be per- 
formed by the dispersive measurement using the non- 
linear Josephson resonato r 12 ' 13 which has the advantages 
of high speed and sensitivity, low backaction, and the ab- 
sence of on-chip dissipation. 

We employed the special conditions g e r r = 7r/2 and 
co e T r = 2ir(p — 1/4) for charge detection. However, the 
general conditions for qubit state readout are given by 



QeJr = 27rg ± 7r/2, 
uJ e T r = lixp ± 7r/2 



(59) 
(60) 



with integers p, q. These conditions result in the rela- 
tion between the coupling constant g e and the Josephson 
coupling energy Ejb as follows, 



9e 



iq±l 
4p± 1 



(61) 



With (4g+l,4p-l) or (4g-l,4p+l) in Eq. ([ST]), we have 
the charge states | j 0o,i( 7r /2)) ~ 1^+1) and |V'2,3(i"/2)) ~ 
|n), whereas with (4g+l, 4p+l) or (4g— l,4p— 1) we have 
|V>o,i(tt/2)) ~ |n) and |^ 2 , 3 (7r/2)) ~ |n + l). In this study 
we use the parameter value, Ejb/Ej=0.1. Hence, the 
coupling strength g e is estimated as g e ~ 10GHz for p = q 
with _Ej=100GHz, which is too strong to be realizable. 
In order to obtain a moderate coupling strength g e we 
need to adopt large p and small q in Eq. (|6"Tj) . We set 
p = 10 and g = 0, and according to Eqs. (|4Tjj) and ((ST)]) 
we estimate the measurement time r r ~ 1 ns, which is 
sufficiently short to maintain qubit coherence during the 
measurement. 



]C-(l/2fr 2 ) cos 2 ruSxr^l with the noise power Sxi- 
Here = tan _1 (A/e) and ?7 c h = ta,n~ 1 (E,jt,/2E c h) for 
flux and charge fluctuation, respectively. Since cos 77a = 
and cos 77^ = at the operating point e = and 
E c h = 0, relaxation is dominant decoherence process 
at this optimal point. The relaxation rate is given by 
1/T-y = X)i(l/2^ 2 ) sin 2 r]iSxi- It is known in experiments 
that Tifi r~j 1/is for flux qubit I 6 and Ti lC h ~ lfxs for trans- 
mon qubilji^, resulting in T\ ~ 0.5/xs. If we increase the 
capacitance of the transmon in the present scheme, the 
relaxation rate can be decreased further. Then, the re- 
laxation time Xi of our qubit can approach that for the 
flux qubit. 

For multi-qubit case the coefficients {oS k \ ft^} may 
be different for different qubit- (k). As shown in Eq. (|57p . 
however, the states \£q(tt/2)) and |£i(7r/2)) produce defi- 
nite charge detection results regardless of the values of a 
and /3. Hence multi-qubit operation and readout of the 
qubit states can also be achieved with these qubits. 

In summary, we propose a readout scheme for the su- 
perconducting flux qubit which is a hybrid of the usual 
three- Josephson junction qubit and the transmon. The 
phase degree of freedom of flux qubit loop and the charge 
degree of freedom of transmon are entangled with each 
other so that the qubit state readout can be achieved by 
detecting the charge number of the transmon. A tt/2- 
rotation of the entangled state by an electric field results 
in the discriminating charge number state which is cor- 
related with the current state of the flux qubit. We show 
that the non-destructive single-shot measurement for the 
flux qubit state can be achieved by detecting the state of 
tranmon. Further, the readout can be performed at an 
optimally biased point with respect to both the magnetic 
field and gate voltage. The fidelity of qubit state readout 
is shown to be sufficiently high. 



IV. DISCUSSIONS AND SUMMARY 
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